Understanding the Unique Structural and Electronic Properties of SrFe02 
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We report a first-principles study of SrFe02, an infinite- layer oxide with Fe atoms in a perfect 
square-planar coordination down to essentially K. Our results reveal this striking behavior relies 
on the double occupation of the 3(i^2 orbitals of high-spin Fe^^. Such an electronic state results from 
the hybridization of iron's Sd^s and 4s orbitals, which enables a large reduction of the intra-atomic 
exchange splitting of the z^ electrons. The generality of the phenomenon is discussed. 



Tsujimoto et ail. [l] have recently reported on SrFe02, 
a novel material that exhibits a wealth of surprising prop- 
erties. Most interestingly, SrFe02 presents Fe^+ ions in a 
square-planar coordination, forming unprecedented infi- 
nite Fe02 layers that remain undistorted down to at least 
4.2 K (see Fig. [T]). In addition, it undergoes a transition 
to an anti- ferromagnetic (AFM) state at Tn=473 K, a 
surprisingly high temperature given the seemingly two- 
dimensional nature of the compound. The physical origin 
of these features, which also appear in similar compounds 
studied more recently [2], remains unknown. Here we re- 
port a first-principles study, based on the "LDA+U" ap- 
proach to Density Functional Theory (DFT), that con- 
firms the peculiar behavior of SrFe02 and reveals its elec- 
tronic origin. 

Basic properties- We begin by showing that our DFT 
methods, described in [s!], reproduce the main results 
reported in Ref. 'jj. To investigate the stability of the 
highly-symmetric phase of SrFe02, we first determined 
the structural parameters that minimize the energy of 
the crystal under the constraints that the P4/mmm space 
group is preserved and the magnetic order is the AFM 
one found experimentally. Our results, shown in Table [H 
are in good agreement with the measurements. We then 
computed the elastic constants and the phonons that are 
compatible with the cell of Fig. [TJ The results, sum- 
marized in Table [H indicate that the P4/mmm phase 
of SrFe02 is stable upon symmetry-lowering structural 
distortions, thus confirming it as the ground state of 
the material at K. It is worth noting that, while ob- 
viously anisotropic, the computed elastic constants and 
phonon frequencies present values that are acceptable for 
a three-dimensional (3D) crystal (for SrFeOs we obtained 
Ci 1=300 GPa and Ci2=140 GPa), which suggests this 
novel infinite- layer phase is robustly (meta-) stable. 

We then studied the magnetic interactions in SrFe02. 
We considered all the spin arrangements compatible with 
the cell in Fig. [H and confirmed that the experimen- 
tally observed AFM pattern (sketched in the figure) is 
the magnetic ground state. Further, the Fe^+ ions are 
in a high-spin configuration with a magnetic moment of 
about 3.6 /iB, in agreement with experiment. We also 



FIG. 1: (color online) 16-atom SrFe02 cell considered in most 
of our calculations. The ground-state magnetic structure is 
sketched; the lattice and exchange constants are defined. 



TABLE I: Computed properties of the P4/mmm phase of 
SrFe02 (measured lattice parameters in parenthesis). 

Structure: a=6=4.01 A (3.99 A); c=3.42 A (3.47 A) 
Elastic constants (GPa): 

B=125, Cii=C22=317, C33=164, Ci2=78, 

Cl3 = C23 = 51, (744 = ^55 = 95, C66 = 125 

Infra-red active modes (cm~^): 

159 and 276 (A2u); 183, 301, and 531 (^u) 
Magnetic couplings (meV): Ji = 1.75, J2=6.58, J3=— 0.26 



computed the exchange parameters of a model Hamilto- 
nian H = ^/'^^ij JijSiSj^ with |5'|=2; the results are 
given in Table [H following the notation defined in Fig.[Tl 
We found that the in-plane interactions are considerably 
stronger than those along the c direction; yet, the mate- 
rial is far from being two-dimensional from the magnetic 
point of view. Indeed, a simple (admittedly rough) mean 
field estimate based on our computed J's renders a Neel 
temperature above 1000 K, and the energy scale for the 
magnetic couplings between Fe02 planes, which we can 
estimate as S'^(2Ji — 8J3), is found to be around 260 K. 



Thus, our results confirm the three-dimensional charac- 
ter of magnetism in SrFe02, and are compatible with 
the relatively high Tn observed. We checked that values 
of the U parameter within the 3-6 eV range, which is 
typical for Fe compounds, do not affect this conclusion. 

Electronic structure- Based on their ^^Fe Mossbauer 
spectroscopy results showing a temperature-independent 
quadrupole splitting of about 1 mm/s, and in accordance 
with the usual picture from crystal-field theory, the au- 
thors of Ref. 1 assumed the {xz^yz)^{xy)^{z'^)^{x'^ — v'^Y 
electronic configuration for the 3(i electrons of Fe^+. 
However, in such conditions a Jahn- Teller distortion is 
expected to lift the degeneracy associated to the (xz^yz) 
orbitals, at odds with the experimental results. Our cal- 
culations readily clarify this question. As shown in Fig. [21 
the DFT electronic structure of SrFe02 does not agree 
with the one proposed in Ref. [1]. We found that the sixth 
3d electron of high-spin Fe^+ goes to the z'^ minority-spin 
orbital. At variance with the (xz^yz) 3d orbitals of Fe, the 
z'^ state is not degenerate with any other and, thus, its oc- 
cupancy cannot cause any symmetry-lowering distortion 
that reduces the energy. Hence, this electronic configu- 
ration is compatible with the stability of the P4/mmm 
phase. Further, we also checked the calculation results 
are compatible with the Mossbauer data. We address 
this point below, after having discussed the electronic 
structure of SrFe02 in more detail. 

Discussion- Our DFT calculations provide us with a 
coherent picture of SrFe02, accounting for the structural 
and magnetic properties of the material. Yet, that pic- 
ture relies on a result, the double occupation of the z'^ 
orbitals of Fe^+, which may seem counterintuitive based 
on the common knowledge in molecular and solid-state 
chemistry. We will thus focus on this key issue: Why are 
the z'^ orbitals of Fe^+ doubly occupied in SrFe02? 

Let us start with a word of caution: Many apparently 
satisfying answers to the above question do not resist 
a rigorous analysis. For example, it may be tempting 
to postulate that the double occupation of the z"^ or- 
bitals is determined by the fact that, as compared with 
the other Sd's, they suffer from negligible Coulomb re- 
pulsion because they find no 0^~ ions along the direc- 
tion perpendicular to the Fe02 plane. This kind of argu- 
ment has been invoked in several occasions [8,, M ^^^ ^^ 
supported by extended Hiickel calculations we have per- 
formed. However, it is severely flawed: It implies that 
the z'^ orbitals should be the lowest-lying ones for both 
the majority- and minority-spin electrons of high-spin 
Fe^+, which is in evident disagreement with the results in 
Fig. [21 Ultimately, the failure of this argument suggests 
that an independent-electron picture is not adequate to 
explain the electronic structure of SrFe02. 

To address the above question, we proceeded by elim- 
inating possibilities in a systematic way. Our key ob- 
servations are the following. As already discussed, our 
calculations show SrFe02 should be regarded as a 3D 
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FIG. 2: (color online) Panel (a): Computed electronic struc- 
ture of SrFe02 around the Fermi energy, which is set to eV. 
The partial Density of States (pDOS) corresponding to the 
3d and 4s orbitals of a particular Fe atom is indicated. The 
dashed line depicts the total DOS divided by a factor of four. 
Panels (b) and (c): Representative F-point eigenstates with 
dominating z'^ character, indicated with arrows in the band 
structure of panel (a). The plotted isosurfaces correspond 
to the indicated values of the squared modulus of the wave 
function. 



crystal, and one can thus conjecture that interactions 
along the c direction, particularly between neighboring 
Fe atoms, might play a role in stabilizing the doubly- 
occupied z'^ orbital. In order to investigate this possi- 
bility, we considered the limit case in which such inter- 
actions are totally suppressed. More precisely, we stud- 
ied an idealized system composed of an isolated Fe02 
layer, fixing the number of electrons in the calculation so 
as to keep the Fe and O atoms in the ionization states 
they have in SrFe02. We denote this virtual 2D crys- 
tal as "Fe^+(0^~)2 +bg", where "bg" refers to the ho- 
mogeneous charge-compensating background introduced 



to stabilize the material. Interestingly, the partial Den- 
sity of States (pDOS) obtained for "Fe2+(02-)2 +bg" is 
essentially identical to the result for SrFe02 in Fig. [21 
indicating that out-of-plane interactions play no role in 
determining the electronic structure of SrFe02. Recog- 
nizing that the significance of the results with a charge 
background may be questionable, we also run calcula- 
tions for SrFe02 with a significantly increased c lattice 
constant and ratified this conclusion. 

Next, we turned our attention to the interactions 
within the Fe02 planes. We ran calculations for slightly 
increased and decreased values of the a=h lattice con- 
stants, and found that the most significant changes in 
the electronic structure pertain to the x^ — y'^ electrons 
of iron. This is not surprising, as these are the elec- 
trons forming cr bonds with the oxygen atoms. In what 
regards the z^ electrons, we observed no effect worth not- 
ing, which clearly indicates that the Fe-0 interactions are 
not the key to understand the double occupation of the z^ 
orbitals. In particular, these results rule out the possibil- 
ity that TT-donor interactions with the oxo ligands, which 
would destabilize the xz and yz orbitals, are at the ori- 
gin of the double occupation of z^. These conclusions are 
ratified by a standard Fe-0 overlap population analysis. 
We also checked that the electronic structure of SrFe02 
is essentially independent of the spin arrangement, and 
we can thus neglect inter-atomic exchange effects in this 
discussion. 

We are thus led to conclude the double occupation of 
the z'^ orbitals must be determined by some sort of intra- 
atomic mechanism associated to the iron atoms. Indeed, 
this is the correct perspective to understand the elec- 
tronic structure of SrFe02. Let us note a well-known 
effect that appears to be critically important in this com- 
pound: As a result of the tetragonal point symmetry at 
the Fe site, the z'^ orbital is fully symmetric and can hy- 
bridize with the formally-empty Fe-45 state. In fact, the 
pDOS of Fig. [2] clearly shows such a hybridization occurs, 
being particularly important for the majority spin. As 
illustrated in panels (b) and (c) of Fig. [21 the majority- 
spin wave functions with a significant As-Sd^^ mixing 
are characterized by (i) a shrinking of the z^-like cen- 
tral ring, and (ii) a deformation, and greater relative 
weight, of the lobes extending out-of-plane. In contrast, 
the minority-spin electrons display a nearly pure z'^ char- 
acter. Note also that the very different structure of the 
z'^ pDOS for the majority (roughly, triple-peaked) and 
minority (single-peaked) spins is consistent with this dif- 
ferent mixing with the 45 orbitals. Of course, this pecu- 
liar behaviour refiects a mechanism by which the crys- 
tal minimizes its energy: Having different spin-up and 
spin-down 2;^-like electrons results in a reduction of their 
Coulomb repulsion and, accordingly, of the corresponding 
intra-atomic exchange spliting A. Indeed, the relatively 
small A associated to the z'^ electrons is obvious from 
Fig. [21 For a given orbital type, A can be estimated as 
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FIG. 3: (color online) Computed energy of various electronic 
solutions of the "Fe^^(0^~)3 +bg" virtual system (see text) 
as a function of the distance d-peO' between the central Fe 
atom and the apical oxygens 0'. 



the difference between the centers of mass of the spin- up 
and spin-down bands; our calculations indicate that A 
is about 2.9 eV for the z'^ orbitals, while the other Sd's 
present A's of more than 6 eV. 

We have thus identified the mechanism that allows 
SrFe02 to reduce its energy by adopting an electronic 
configuration with doubly-occupied z'^ orbitals. Inter- 
estingly, for this mechanism to be active we only need 
to have the appropriate local symmetry at the Fe site, 
which implies that the newly discovered effect might be 
quite general. To check this, we performed the following 
computational experiment: We considered the 3D crystal 
"Fe^+(0^~)3 +bg" sketched in Fig.[3l fixing the number 
of electrons in order to retain the Fe-3d^ configuration. 
For simplicity, we assumed a ferromagnetic order, which 
we checked has no effect on the electronic configuration of 
Fe. We fixed the in-plane lattice constant at the SrFe02 
value, and monitored the evolution of different electronic 
solutions [lO[ as the cubic crystal ((iFeO'=2 A) under- 
goes a tetragonal distortion. As shown in Fig. [3l only 
one electronic solution is found in the cubic case: the 
minority-spin 3d electron occupies a combination of t2g 
orbitals (from which the xy orbital is selected as soon 
as the tetragonal distortion occurs), closing the gap of 
the crystal. In contrast, for the largest values of d^eO' 
considered, several solutions occur, the ground state cor- 
responding to doubly-occupied z'^ orbitals. Interestingly, 
the solution with doubly-occupied z'^ already exists at 
<^FeO'=2.14 A; yet, the "x?/-to-z^" crossover only occurs 
at (iFeO'=2.50 A. These results show a mere symmetry 
breaking will not cause the double occupation of the z'^ 
orbitals; rather, a relatively large tetragonal distortion, 
which results in both a more anisotropic environment for 
the Fe atom and a reduction of the Coulomb repulsion 
for the z'^ orbitals, is required. Finally, we also checked 



that changes in the in-plane lattice constant have a minor 
effect in the xy-to-z'^ crossover. 

It is not easy to estimate the universality of this mech- 
anism for the double occupation of the z'^ orbitals. It 
may seem surprising, for example, that such an elec- 
tronic solution is not found in molecules with Fe^+ cen- 
ters. However, after reviewing the literature for some 
well-known cases [11], we have the impression that the 
electronic structure of iron is often conditioned by the 
nature of the ligands or by a reduced atomic symme- 
try that precedes any electronic effect associated to the 
metal atom itself. This suggests such systems may not 
be appropriate to look for effects like the one discussed 
here. As for crystalline materials, we only know of one 
candidate to present this kind of effect: the mineral gille- 
spite (BaFeSi40io), which contains Fe^+ ions in a square- 
planar coordination and displaying a high-spin configura- 
tion with doubly-occupied z'^ orbitals [sl]. An appealing 
possibility, though, would be to synthesize other com- 
pounds of the SrFe02 family. Indeed, we found that 
crystals like RbCo02 and YFe02 would present a similar 
pattern of electronic levels with doubly-occupied z'^ or- 
bitals. More precisely, our calculations predict RbCo02 
displays a Co-Sd^ high-spin configuration with a doubly- 
occupied z^; thus, this material might remain in the 
high-symmetry P4/mmm phase down to K, exactly as 
SrFe02. On the other hand, YFe02 exhibits a Fe-3d^ 
high-spin configuration where, in addition to having a 
doubly-occupied z^, the last 3d electron goes to a (xz^yz) 
state that closes the band gap; hence, this crystal should 
undergo a Jahn- Teller distortion. 

Finally, let us discuss the compatibility between the 
value of the quadrupole splitting (QS) obtained from 
^^Fe Mossbauer spectroscopy, which is about 1 mm/s, 
and the double-occupation of the z'^ orbitals, which is 
expected to result in QS values of 2-3 mm/s. We should 
begin by noting that SrFe02 is not the only compound 
exhibiting this seemingly contradictory behaviour. In- 
deed, the Fe^+ atoms in the above mentioned gillespite 
have been reported to display QS values below 0.6 mm/s 
in the 80-650 K temperature range [12]. As for SrFe02, 
we computed the electric field gradient at the Fe site 
and the corresponding QS, obtaining a value of about 
0.8 mm/s. We also checked that varying U between 
and 6 eV renders QS results in the 0.6-0.9 mm/s range, 
the greater values corresponding to the smaller [/'s. This 
clearly shows that having an electronic structure with 
doubly occupied z'^ orbitals is not incompatible with the 
relatively small QS values measured. As to the mech- 
anism that causes such a reduction of the observed QS 
with respect to its expected value, two possibilities can 
be mentioned: (i) that the lattice contribution to the 
QS is comparable to the one coming from the 3d elec- 
trons of iron [12], which would lead to a cancelation as 
their signs are known to be opposite, or (ii) that the 
splittings associated to iron's 3d electrons are unusual in 



SrFe02. Interestingly, our results seem to support the 
latter possibility: The spread in energy of the majority- 
spin z'^ and x^ — y'^ electrons (see Fig. [2^) will imply a 
relatively large spacial spread of the corresponding Wan- 
nier functions [l3| , which might result in unusually small 
quadrupole splittings. However, a detailed analysis of 
this issue is not trivial and goes beyond the scope of this 
paper. 

Summary- We have used first-principles methods to 
show that SrFe02, an unusual high-spin Fe^+ layered 
phase recently reported, is free from structural insta- 
bilities and exhibits strong magnetic interactions. We 
find that the hybridization of iron's 3dz'2 and 45 orbitals, 
which results in a large reduction of the intra-atomic ex- 
change splitting associated to the z'^ electrons, is the key 
feature to understand the structural stability of the ma- 
terial. 
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